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摘要 : 望远镜 自动 化 观测 编排 ,是 望远镜 观测 控制 系统 实现 自主 观测 必 不 可 少 的 组 成 部 
分 。 云 南天 文人 台 40 米 射 电 望 远 镜 需要 对 观测 目标 制定 中 长 期 和 短期 的 纲要 ， 目 前 对 短期 的 
目标 进行 纲要 的 编制 时 ， 主 要 依赖 观测 人 员 根 据 经 验 手动 来 实现 ， 这 种 方式 不 仅 低 效 耗 时 ， 
更 有 可 能 会 导致 结果 显著 地 偏离 最 优 观 测 方案 , 不 能 充分 利用 望远镜 宝贵 的 观测 时 间 。 本 文 
使 用 混合 整数 线性 规划 模型 ， 分 析 了 云南 天 文 台 40 米 射 电 望 远 镜 的 短期 观测 目标 的 问题 ， 
定义 了 短期 观测 目标 调度 中 的 限制 条 件 ， 以 最 小 化 望远镜 的 寻 址 时 间作 为 目标 函数 ,使 用 数 
学 优化 来 解决 射电 望远镜 短期 观测 目标 编排 问题 ， 并 在 云南 天 文 台 40 米 射 电 望 远 镜 上 进行 
了 测试 ， 能 较 好 地 满足 目前 短期 观测 目标 编排 的 需求 。 
关键 词 : 观测 编排 ; 线性 规划 ; 射电 望远镜 ; 自主 观测 
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1 引 


为 了 减少 电磁 干扰 对 观测 的 影响 ， 射 电 望 远 镜 台 站 都 选择 在 高 海拔 ， 人 迹 稀少 的 地 区 ， 
但 这 同时 也 限制 了 人 工 现场 干预 观测 能 力 。 同 时 当前 望远镜 正 朝 着 巨型 化 、 复 杂 化 发 展 ， 望 
远 镜 的 正常 运转 需要 多 种 组 件 协同 工作 叫 , 依赖 操作 人 员 手 工控 制 已 变 得 越 来 越 困难 。 因此 ， 
自主 观测 已 是 每 个 望远镜 观测 控制 系统 的 共同 需求 , 自主 观测 是 指 根据 观测 计划 自动 进行 调 
度 的 观测 模式 ， 而 观测 计划 的 编排 是 望远镜 实现 自动 化 观测 控制 的 前 提 。 
观测 计划 编排 已 被 广泛 研究 , 但 主要 集中 在 长 期 观测 申请 调度 上 , 通常 是 望远镜 科学 委 
员 会 根据 科学 目标 ,赋予 不 同 的 观测 申请 不 同 的 科学 价值 权重 后 , 通过 优化 方法 给 出 长 期 的 
计划 编排 证。 而 对 于 云南 40 米 射电 望远镜 一 个 重要 需求 是 ， 给 出 一 系列 的 观测 目标 ， 使 望 
远 镜 尽 可 能 在 短 时 间 内 完成 所 有 观测 目标 的 观测 ,通常 是 一 天 内 完成 几 十 颗 源 的 短期 观测 编 
HE, 显然 目前 在 其 他 望远镜 系统 的 长 期 观测 编排 方法 并 不 能 满足 云南 40 米 射 电 望 远 镜 短 区 
编排 需求 。 

不 同 于 长 期 观测 编排 , 短期 编排 需要 给 出 观测 目标 具体 的 开始 和 结束 观测 时 间 , 而 时 间 
变量 是 连续 的 ， 因 此 短期 编排 属于 连续 最 优化 问题 。 观 测 目 标的 位 置 随 着 时 间 不 停 变 化 ， 而 
不 同 的 观测 目标 一 天 内 在 望远镜 视 场 中 出 现 的 时 间 有 限 , 同时 还 要 考虑 从 一 个 观测 目标 转动 
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到 另外 一 个 目标 , 由 于 望远镜 本 身 机 械 问 题 , 可 能 出 现 的 限 位 以 及 观测 目标 能 否 获得 较 好 的 
数据 质量 , 这 些 都 是 短期 编排 区 别 于 长 期 编排 的 地 方 ， 也 是 优化 的 难点 。 本 文 首先 回顾 了 当 
前 主要 的 观测 编排 方法 , 接着 介绍 将 望远镜 的 源 的 编排 转换 成 混合 整数 线性 规划 问题 中 的 决 
策 变 量 、 目标 函数 和 约束 条 件 的 定义 , 给 出 了 该 方法 在 云南 40 米 射 电 望 远 镜 上 的 测试 结果 ， 
最 后 是 本 文 工 作 的 总 结 ， 并 讨论 了 云南 40 米 射电 望远镜 观测 编排 未 来 的 研究 工作 。 


2 相关 工作 


调度 问题 已 在 电力 、 运 输 、 航 空 、 计 算 机 以 及 工业 生产 等 领域 得 到 了 广泛 的 研究 ， 望 远 
镜 观测 目标 的 编排 正 是 该 问题 在 天 文 领域 的 实际 应 用 。 哈 勃 空 间 望 远 镜 的 Spike 系统 自 1990 
年 起 投入 使 用 ， 提 出 的 核心 算法 为 迭代 修正 搜索 算法 外 。 迭 代 修正 搜索 算法 的 核心 思路 是 将 
望远镜 时 间 分 配 问题 抽象 成 一 个 约束 优化 问题 ,随后 通过 三 个 步骤 进行 求解 : 首先 通过 某 种 
启发 式 算法 构造 初始 解 ,该 初始 解 可 以 违反 约束 条 件 ， 其 使 用 的 启发 式 算法 可 以 根据 具体 需 
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N 求 进行 选择 ,但 启发 式 过 程 均 基于 对 约束 程度 的 量化 并 且 选 择 容易 计算 的 量化 方式 ,随后 再 
= 利用 启发 式 的 修正 算法 来 减少 约束 的 违反 。 最 终 通过 从 时 间 分 配 表 中 去 除 存在 冲突 的 活动 来 
© MERU, 


绿 岸 射电 望远镜 (Green Bank Telescope, GBT) 采用 了 NRAO 开发 的 一 个 自动 动态 调 
度 系统 名 ， 其 核心 思想 是 基于 一 种 排序 的 算法 来 确定 观测 计划 之 间 的 优先 次 序 。 文 献 描述 了 
GBT 使 用 的 包括 预报 天 气 条 件 、 观 测 效率 、 目 标 坐 标 等 变量 的 排序 算法 ， 以 及 调度 观测 中 
使 用 的 一 些 特定 变量 。 

撒 丁 岛 射 电 望 远 镜 (Sardinia Radio Telescope, SRT) 中 参考 了 GBT 的 动态 调度 系统 ， 

考虑 了 长 期 和 短期 编排 问题 ， 在 其 模拟 结果 中 ， 以 48 小 时 为 周期 ， 每 一 个 小 时 作为 一 个 时 
间 槽 。 每 个 观测 计划 被 赋予 一 个 科学 价值 ， 占 用 2 到 6 小 时 ， 使 用 模拟 退火 算法 优化 观测 编 
排 ， 加 入 了 水 汽 和 天 气 条 件 等 限制 条 件 。 
500 米 口径 球面 射电 望远镜 (Five-hundred-meter Aperture Spherical radio Telescope, FAST) 
中 以 最 小 化 寻 址 代价 和 最 大 化 科学 价值 作为 目标 函数 ， 限 制 条 件 考 虑 了 源 的 可 见 性 ， 射 频 干 
扰 等 因素 , 同时 使 用 了 边缘 检测 的 方式 减少 问题 空间 ， 加 快 优 化 速度 。 使 用 了 中 国 虚拟 天 文 
台 的 在 线 观 测 提案 系统 ， 本 质 上 还 是 中 长 期 的 观测 编排 。 
其 他 的 望远镜 系统 的 观测 编排 也 大 都 是 观测 提案 的 编排 ， 比 如 基于 进化 算法 的 LOFAR 
观测 调度 如， 基于 模拟 退火 和 动态 规划 实现 的 ALMA 观测 调度 外 。 而 云南 40 米 射 电 望 远 镜 
作为 单 碟 射 电 望 远 镜 ， 在 短期 的 观测 编排 中 ， 需 要 在 更 短 的 时 间 模 内 ， 按 目标 源 进行 分 配 ， 
而 不 是 分 配 观测 提案 。 而 现 有 的 中 长 期 的 观测 提案 的 编排 方法 ， 显 然 无 法 解决 按 目标 源 进行 
观测 时 间 分 配 的 需求 。 
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3 模型 设计 


3.1 问题 描述 


云南 40 米 射 电 望 远 镜 短期 观测 目标 编排 需求 通常 是 在 一 天 《或 几 个 小 时 ) 内 ， 给 定 一 
系列 的 观测 目标 信息 ， 包 含 观测 目标 的 名 称 (J2000 纪元 )、 位 置 和 所 需 观测 时 长 ， 如 表 1 


所 示 。 
#1. 云南 40 米 射电 望远镜 观测 源 信息 样 例 

观测 源 名 称 IRA CRA) 赤 纬 (DEC) 观测 时 长 分钟》 

0034-0721 00 34 08.8703 -07 21 53. 409 48 

0139+5814 01 39 19. 7401 +58 14 31.819 48 

1559-4438 15 59 41. 5261 -44 38 45. 9018 24 

1745-2900 17 45 40. 190005 -29 00 30. 37032 48 

短期 编排 的 结果 是 得 出 每 个 目标 开始 和 结束 的 观测 时 间 ， 并 使 望远镜 的 利用 率 最 大 化 。 
同时 需要 考虑 : 


1) 所 有 观测 目标 都 被 安排 观测 。 每 颗 目 标 源 的 观测 时 长 要 求 集中 在 24 到 48 分 钟 ， 需 
保证 编排 时 间 满 足 观测 时 长 要 求 。 

2) 分 配 时 间 粒 度 。 时 间 分 配 粒度 太 大 ,会 造成 时 间 片 的 浪费 。 如 果 时 间 分 配 粒 度 太 小 ， 
会 造成 参数 空间 的 膨胀 ， 直 接 影响 最 优化 求解 的 性 能 

3) 望远镜 硬件 的 限 位 条 件 。 在 观测 过 程 中 源 的 位 置 不 断 变 化 ， 而 望远镜 的 俯仰 角 只 能 
在 一 定 的 范围 转动 ， 因 此 须 满足 在 源 的 编排 时 间 范 围 内 ,望远镜 都 可 以 转动 到 源 目标 , 也 即 
是 满足 望远镜 硬件 的 限 位 条 件 

4) 最 优 观 测 时 间 。 根 据 40 米 望 远 镜 的 历史 观测 经 验 ， 天 线 俯仰 在 40-60 度 内 进行 观测 
比较 好 ， 这 时 候 增益 较 大 ， 且 增益 变化 不 大 ， 观 测 源 受 地 面 噪声 相对 较 小 ， 观 测 获 取 的 数据 
质量 也 就 更 好 。 因此 编排 时 , 优先 考虑 这 个 位 置 范 围 的 源 , 源 的 位 置 高 度 可 用 Airmass 衡量 ， 
Airmass 值 可 近似 于 1/sin a, 其 中 a 为 源 的 高 度 角 。 如 图 1 所 示 , 显示 了 2020 年 12 月 21 H, 
源 2055+3630，0837-4135，0358+5413 的 Airmass 随时 间 变 化 的 曲线 。 
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Figure 1. The Airmass of 2055+3630, 0837-4135, 0358+5413 at different times on December 21, 2020 

根据 云南 40 米 射电 望远镜 短期 观测 目标 编排 需求 和 限制 条 件 ， 假 设 总 共有 天 个 观测 目 
标 , 一 个 观测 目标 在 观测 过 程 中 需要 考虑 目标 的 观测 时 长 以 及 望远镜 从 第 一 颗 源 的 跟踪 状态 
切换 到 第 二 颗 源 的 跟踪 状态 所 消耗 的 时 间 〈Slew Time )， 因 此 ， 一 个 观测 目标 上 占用 的 时 间 
窗 为 转向 时 间 Ts+ 观 测 时 间 T。， 可 表示 为 ; 


Tear +h a) 
其 中 转向 时 间 依 赖 上 个 观测 目标 结束 时 望远镜 的 位 置 〈 方 位 、 俯 仰角 度 )。 如 果 是 首 个 
观测 目标 ， 则 望远镜 的 转向 时 长 根据 望远镜 初始 的 方位 角 和 高 度 角 进行 计算 。 
3.2 决策 变量 
在 望远镜 短期 观测 编排 中 ,时 间 是 最 主要 的 变量 , 时 间 直 接 决 定 了 源 的 位 置 ， 也 间接 影 
响 源 的 观测 时 长 。 在 云南 40 米 的 短期 目标 编排 中 ， 目 标 源 观测 时 长 的 单位 是 分 钟 ， 我 们 将 
时 间 的 分 配 粒度 设置 为 分 钟 ， 整个 观测 计划 的 总 时 间 M 为 24 x 60 x d 分 钟 , 其 中 d 为 编排 
计划 的 天 数 。 
首先 对 于 每 个 观测 目标 {Tj, jE K}， 定 义 目标 7) 的 开始 观测 时 间 为 台 € [0,M]， 结 束 观 


测 时 间 为 € [0,M]， 在 优化 实现 中 ， 为 加 快 求解 速度 ,如 和 ty 范围 取 目 标 7 第 一 个 可 用 时 间 


槽 的 开始 时 间 和 最 后 一 个 可 用 时 间 槽 的 结束 时 间 。 我 们 取 Airmass 的 范围 为 [2, 10]， 可 计算 


获得 每 个 观测 目标 在 观测 总 时 间 内 的 s 个 可 观测 的 时 间 模 ,如 图 2 所 示 , 每 个 观测 源 在 一 天 
内 通常 有 1-4 个 可 观测 的 时 间 槽 ， 单 个 目标 在 所 有 可 观测 时 间 覃 的 Airmass 均 在 [2, 10]. xe 


义 二 进 制 变量 Ya = {5， 表 示 目 标 Tj 是 否 在 其 第 nn 个 可 观测 的 时 间 村 内 被 调度 。 
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图 2. JH 0922+0638 和 源 1136+1551 在 2020 4 12 H 21 AAAS AA Ta 


Figure 2. The available time slots of 0922+0638 and 1136+1551 on December 21, 2020 


观测 目标 编排 重要 的 点 就 是 编排 源 的 观测 顺序 ， 定 义 二 进 制 变量 yi fo 表示 望 远 


镜 是 否 完成 源 7, 的 观测 后 转动 到 源 7 o 
定义 目标 源 的 距离 矩阵 为 D， 其 中 目标 Ti(@1, 61) 和 目标 Tj(qz, 62) 间 的 距离 di， 使 用 大 
ILS 
dj; = arccos [sin 5; sind; + cos 61 cos 62 cos (a — a2)] (2) 
云南 40 米 望 远 镜 根据 经 验 值 ， 最 大 的 源 目 标 切 换 时 间 约 为 6 ah, AUR BRIN Ta 
阵 近 似 于 : 
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Sx6x D (3) 


其 中 万 为 距离 矩阵 了 的 0 均值 归 一 化 矩阵 。 


3.3 目标 函数 


基于 最 大 化 望远镜 利用 率 ， 尺 量 减少 总 的 望远镜 源 切换 时 间 的 原则 ， 目 标 函 数 定义 为 : 


min dies dies Vij G = K) + YijSij (4) 


3.4 约束 条 件 


根据 云南 40 米 射 电 望 远 镜 短期 观测 目标 编排 的 问题 描述 ， 在 最 优化 目标 函数 的 过 程 中 
要 满足 以 下 约束 : 
约束 一 : 满足 每 个 观测 目标 的 观测 时 长 的 要 求 。 即 满足 结束 调度 时 间 减 去 开始 调度 大 于 
于 观测 时 长 Pj。 


s 


e _ +8 : ; 
t-t; 2 P; |jeKk (5) 


约束 二 : 每 个 目标 有 不 止 一 个 可 观测 时 间 槽 , 需要 保证 开始 时 间 和 结束 时 间 落 在 同一 个 
时 间 模 内， 如 图 3 所 示 ， 观 测 源 0922+0638 出 现 的 编排 1 并 不 满足 约束 二 ， 结 束 时 间 并 未 
落 在 时 间 槽 内。 编排 2 满足 约束 二 ， 开 始 时 间 和 结束 时 间 落 在 同一 个 时 间 槽 内 。 


00:00 ”092240638 01:54 15:36 09220638 19:23 21:32 092210638 24:00 
3， 约 束 二 条 件 限制 示意 图 


Figure 3. The diagram of condition restriction under the second constraint 


也 即 是 如 果 目 标 帮 的 开始 观测 时 间 为 她 大 于 等 于 一 个 时 间 档 的 开始 时 间 ， 则 结束 时 间 妃 


应 小 于 等 于 这 个 时 间 模 的 结束 时 间 。 


站 (6) 


其 中 N 是 目标 Tj 的 可 观测 时 间 模 的 数量 , LMS n SSS TPS Ti], Lem n 个 时 间 


槽 的 结束 时 间 。4 = {09， 表 示 目标 7 是 否 在 其 第 nn 个 可 观测 的 时 间 村 内 被 调度 ， 且 需 满 足 


olin = 1， 表 示 目 标 7) 必 须 在 其 中 一 个 时 间 槽 内 被 调度 。 
约束 三 ， 保 证 每 个 观测 目标 开始 观测 时 间 和 观测 结束 时 间 没 有 相互 重合 。 


if yy=1 > t > tt Sy (7) 


— = 1, 表示 完成 源 T; 的 观测 后 转动 到 源 Tj, m BEN EAT HSE GY TA AP PRT 
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4 实验 


实验 数据 来 自 云南 40 米 射 电 望 远 镜 的 真实 观测 源 ， 短 期 编排 的 时 间 区 间 为 1 天 ， 共 有 
观测 源 24 颗 。 每 颗 源 所 需 观 测 时 长 在 24-48 分 钟 。 实 验 的 硬件 环境 为 1 台 配 置 为 Intel Xeon 
CPU E5-2660 v4 CPU @ 3.4GHz 处 理 器 (56 核 )，512GB RAM 的 高 性 能 服务 器 。 求 解 器 使 用 
Gurobi 9.1， 以 2020 年 12 年 21 日 为 例 ， 如 图 4 所 示 ，24 颗 目 标 源 被 全 部 编排 ， 并 不 存在 
时 间 交 叉 的 情况 ， 经 过 验证 也 都 在 其 可 观测 的 时 间 槽 内 。 


O ROY 


O O O O 
0° se abo 


Minutes since 2020-12-21 00:00:00 
Al 4. 24 颗 观测 源 在 2020 年 12 月 21 日 的 观测 时 间 分 配 


Figure 4. The observation time scheduling results of the 24 sources on December 21, 2020 
Gap 即 百分比 界 差 , 是 求解 大 规模 混合 整数 规划 问题 时 非常 重要 的 参数 ， 可 视 为 在 去 除 
问题 参数 整数 取 值 的 限制 后 ， 问 题 松 弛 解 与 最 优 解 之 间 的 绝对 差距 : 


| 
LP 


Gap = | (8) 


OPT 为 当前 最 优 解 , LP 为 去 掉 整 数 约束 后 的 松弛 最 优 解 。 实验 环境 中 , 使 用 56 个 线程 ， 
模型 可 在 300 秒 内 将 Gap 控制 在 0.1 以 下 。 


5 总 结 


本 文 利用 线性 规划 方法 对 云南 天 文 台 40 米 射 电 望 远 镜 的 短期 观测 目标 编排 问题 进行 了 
建 模 , 详 述 了 在 短期 观测 编排 过 程 中 过 到 的 限制 条 件 , 利用 非 线 性 混合 整形 线性 规划 解决 了 
这 一 约束 优化 问题 。 实 验 结果 显示 了 该 方法 的 可 用 性 。 

短期 观测 目标 编排 问题 中 ， 机 会 源 的 临时 插入 观测 ， 整 个 观测 计划 就 需要 重新 编排 ， 因 
此 求解 速度 也 是 短期 观测 目标 编排 中 的 重要 考量 。 未 来 研究 使 用 其 他 的 优化 算法 进一步 提高 
编排 性 能 ， 减 少 望 远 镜 等 待 编排 的 时 间 。 同 时 考虑 加 入 目标 源 的 权重 ,在 望远镜 空闲 时 间 对 


权重 高 的 目标 多 次 观测 。 此 外 研究 如 何 利用 已 观测 的 脉冲 星 数据 , 反 演出 天 线 干 扰 情 况 与 方 
位 以 及 俯仰 之 间 的 关系 图 , 在 编排 目标 源 时 优先 干扰 较 小 的 源 。 在 目前 的 模型 中 源 切 换 的 时 
间 使 用 球面 距离 近似 估计 , 不 完全 准确 , 未 来 研究 如 何在 优化 过 程 中 根据 源 的 方位 角 和 高 度 
以 及 望远镜 的 实时 指向 位 置 动 态 计 算 源 的 切换 时 间 。 
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A Short-term Scheduling Method Based on Linear 


Programming for Yunnan 40-meter Radio Telescope 
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2, Yunnan Observatories, Chinese Academy of Sciences, Kunming 650216, China; 
3, Key Laboratory of Applications of Computer Technology of the Yunnan Province, Kunming 650500, China) 

Abstract: The observation scheduling is an indispensable part of the telescope observation 
control system to realize autonomous observation. The 40-meter radio telescope of the Yunnan 
Observatories combines the requirements of long-term and short-term observation scheduling. In 
the current practice of short-term observation scheduling, it mainly relies on the observers to 
manually implement them based on experiences. This method is not only inefficient and 
time-consuming, but also cause the results to deviate significantly from the optimal observation 
plan, and the valuable observation time of the telescope cannot be fully utilized. In this paper, the 
mixed integer linear programming model is used to analyze the problem of the short-term 
observation target of the 40m radio telescope. We define the constraints in the short-term 
observation target scheduling, adapt the minimization of total slew times as the objective 
function, and use mathematical optimization to solve the short-term observation target scheduling 
problem. The proposed method has also been tested on the 40m radio telescope, which can better 

meet the requirements of current short-term observation target scheduling. 
Key words: Observation scheduling; Linear programming; Radio telescope; Autonomous 


Observation 


